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Abstract 

The structure functions bi2{x) of the deuteron are studied within covariant ap- 
proach. It is shown that usual nonrelativistic convolution model result in incorrect 
behavior of this structure functions at small x and violates the exact sum rules. Re- 
alistic calculations are carried out with the Bethe-Salpeter amplitude of the deuteron 
and compared with the nonrelativistic results. 



1 Introduction 

The study of the deep inelastic lepton scattering with polarized targets and beams provides 
refined information about the quark content of hadrons. These days much attention is 
attracted to the nucleon's spin-dependent Structure Functions (SF), g^2 (see e.g. review 
and references therein). The SFs, (71 2, is the simplest example of the spin-dependent SFs, 
which exist for all targets with a nonzero spin, starting s = 1/2. At the same time, hadrons 
with a spin higher than 1/2 have additional spin-dependent SFs [@]. 

The spin-dependent SFs 6i,2(a;)0of spin-1 hadrons have been studied on a few occasions, 
including the vector mesons and deuteron [0, ^, ^ ^ . For mesons some qualitative estimates 
have been done, but "a real understanding of hi{x) at the quark level is not yet available" [Q. 
The only exact sum rules for mesons have been proposed by Efremov and Teryaev IQ: 

1 

Jh{x)dx = 0. (1) 


1 

J h2{x)dx = 0. (2) 


It would be unrealistic to suggest that these sum rules will be experimentally verified for 
the mesons any time soon. However, they are independent of the target, i.e. it is supposed 
to be valid for the spin-1 nuclei as well. The deuteron is the most probable candidate 
for its SFs 61 2 to be measured. Indeed, a number of deep inelastic experiments on the 
polarized deuterons are run or being prepared at the world's best facilities, such as SLAC, 
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^The notation 61,2 is used, following rcf. 
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CERN, CEBAF and DESY. These experiments are usually aiming to extract a neutron SF 
gi, however, in principle, the same data can be used to obtain SFs 61,2- The problem is 
caused only by the yet insufficient accuracy of measurements. Realistic calculations with 
the deuteron wave functions in the Bonn and Paris potentials have been done in ref. 0. 
The preliminary results within the Bethe-Salpeter formalism have also been presented [0. 
It was shown that relativistic calculations differ from nonrelativistic ones, but the reasons 
have not been discussed. The sum rules have not been analyzed in these papers but 

the displayed shape of the SFs makes questionable the validity of the sum rules. 

This letter presents a study of the deuteron SFs &f2(^) within a covariant approach, 
based on the relativistic convolution formalism for the deep inelastic scattering ^, 0, |ll| 
and the Bethe-Salpeter formalism for the deuteron bound state [^, |TB|, T^\. The issue of the 



sum rules is specially addressed. It is shown that the nonrelativistic convolution model is 
not relevant to calculate SFs (^) small x, because this model inevitably breaks at least 
one of the sum rules (|l]) or (Q) . At the same time, the relativistic calculations based on the 
Bethe-Salpeter amplitude for the deuteron are consistent with both of the sum rules. The 
SFs, 

6^2 (•^)' calculated within both relativistic and nonrelativistic approaches. 

2 Structure functions in the relativistic impulse ap- 
proximation 

The nucleon contribution to the deep inelastic scattering of electrons off the deuteron is 
defined by the triangle diagram (see Fig. |I]) |T^, [TT]]: 

{W)m = I J (^Tr {^M{po,p)W{q,pi)^M{po,p){p2 - m)} , (3) 

where pi^2 = Pd/'^^P = (^^d/2±po, Pd is the deuteron momentum, Md is the deuteron 
mass, p = {po,p) is the relative momentum of nucleons, W is an operator appropriate for 



the process (see e.g. and references therein [0, |T5|, 0) ^m(Po?p) is the Bethe-Salpeter 



amplitude for the deuteron with M being the deuteron's total momentum projection. The 
original convolution model for deep inelastic scattering is reproduced by choosing operator 
W in the form: 

^ = ^w;:mpi). (4) 

where m is the nucleon mass, q = {u, 0, 0, —\/v^ + Q^) is the momentum transfer, = —q 
and W^uip, q)^ is the hadron tensor of nucleon: 
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w!!M.p)= (5) 

Ff(a:,g2) 



^-f^Fnx,Q') + [p,~q, 




pq 
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where x = Q'^/{2pq) anf Ff^ are the nucleoli SFs. The small effects of the off-mass-shell 



deformation of the nucleon tensor [jT5|, |T6|, |T^ are not considered in the present letter, because 
these effects do not affect sum rules (|l|)-(@) and do not noticeably change the absolute values 
of the SFs. That is why SFs F/^ in (|^) do not depend on p"^, but and pq. 

One can extract scalar SFs from the hadron of the deuteron. For instance, for F2 one 
gets: 



2MdJ (27r)4 ' \pio + pn 

Tr{^'A,/(po,p)(7o + 73)*m(po,p)(P2 - m)} , (6) 

where xn = (5^/(2mz/) is the Bjorken scaling variable^, i.e. this is x for the on-mass-shell 
nucleon at rest, pio and pis are the time and 3-rd components of the 1-st nucleon momentum. 
Formula (j^) has not been averaged over the projection of the deuteron total momentum, M, 
since in the present form it gives an understanding what the SF b2 is. Indeed, eq. (^) 
gives two independent " SFs", with M = ±1 and M = 0, which are related to the usual 
spin-independent SF, F^, and a new SF, 6^: 

Fi'{x^,Q') = \ E F2''{xr,,Q\M), (7) 

-J M=0,±1 

b2{xN, Q^) = F^{x, Q\ M = +1) - F^{x, Q\ M = 0), (8) 

(xjv, Q^ M = +1) = (x^, Q\ M = -1). (9) 

Of course, instead of (0) and (||) any other linearly independent combination of the functions 
F2^(x, Q^, M) can be chosen. 

Note, SFs F^{x, Q^, M) are independent of the lepton polarization, therefore, both SFs, 
F2 and 6^, can be measured in experiments with an unpolarized lepton beam and polarized 
deuteron target. In view of eq. (P), only one of the SFs F^{x, Q^, M) is needed, in addition 
to the spin-independent F^{x, Q^), in order to obtain b2{x, Q"^). The other SF, 6f , is related 
to the deuteron SF F/^, the same way as 6^ is related to F2 , viz. via eqs. (^, and 6f = 2a;6f . 



3 Singularities of the triangle diagram and sum rules 

It has been previously shown [|, |, |Tl| how a singular structure of the triangle graph 
(Fig. |I]) rules the behavior of the spin-independent SF F^. In particular, it has been found 
that the relativistic impulse approximation satisfies the unitarity and provides the correct 
kinematical region of the variable xn- However, both these properties of the exact covariant 
amplitude are broken in practical calculations, when nonrelativistic wave functions of the 
deuteron are used. In this case one can refer to the argument that such deviations are small, 
and are not important for phenomenology. At the same time, a realistic Bethe-Salpeter 

^Note that the "native" deuteron variable is xd ~ {'m/M£i)xN , however xn is used more often. 
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amplitude of the deuteron serves ideally for a consistent phenomenological application of the 
covariant theory of the processes on the bound nucleons. 

In order to calculate SFs, (H)-® and analyze the sum rules, the singularities of the 
triangle diagram should be explicitly taken into account. To do that, eq. is rewritten as: 

I f d^p N ( 2^ 1 

- — —5 Q 



Tr {(f)M{Po,p){Pi + m){-fo + 73) (pi + m)(f)M{po,p){P2 + m)} , (10) 

where (f>M{Po,p) = {Pi — ™)'^ m{Po,p){P2 — is the Bethe-Salpeter vertex functions of the 
deuteron. 

Analysis of singularities in the complex p2+-plane allows for one analytical integration in 
( PUD \ W^ . Being translated into variables which are used in the present paper, this integration 
is equivalent to picking the residue in the second nucleon pole, P20 = 00 = \/ vn? + or 
Po = Mj:)/2 — uj, in the complex po-plane when both of the following conditions are satisfied: 

Q <uj -p:i< Md. (11) 

The contribution of the region of p beyond (pH]), into integral ([T0|), is zero, i.e. different poles 
cancel each other. Note, that in the required pole pio = M£, — uj. Calculating residue in 
(p!OD, one gets: 

F-(.., g^ M) = -i- / f TT^^^^, Q'] (12) 

Tr{0A/(Po,p)(Pi + "^)(7o + 73)(Pi + "^)0m(Po,p)(P2 + "^)} , 

where the 9-function guaranties the right of conditions (pUf) , the left condition is always 
satisfied. 

It is useful to rewrite (IT^) in the convolution form: 



y 

u 

where "the effective distribution" of nucleons in the deuteron is defined by 



F2^(x^,Q^M)= / dyF^r^M'\fJ\v). (13) 



fT{y) (14) 

1 ^ ''^-Ay-^^^^^\^^yy 



2MdJ (27r)3 V m J ' 2ujMI{Md - 2ujY 

Tr |0m(Po,p)(Pi + "^)(7o + lz){Pi + "^)0m(Po, p)(P2 + m)| 
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Two sum rules can be written down for the effective distribution fM^{y)- 

M]j/m 

fZ"'{y)dy={D\Q\D)M = l. (15) 





Mb jm 



J yf^/''{y)dy={D\{QMr,\D)M = l-SM, (16) 



where Q oc ip{x)'yo4'{x) is the vector charge and {Qn)1^ oc i%l){x)'~^^d^ip{x) is the trace of 
energy-momentum tensor. Eq. ( [T3| ) presents the vector charge conservation generahzed for 
the deuteron states with different M . In spite of such clear physical interpretation, some 
time ago it was a subject of some controversy D|, |TD[. Indeed, the derivation of sum rule 



( p!5|) contains some subtle points and equivalence between it and the expression for the charge 

{D\Q\D)m = ^J ^,li:Tr{^M{Po,Pho^M{Po,p){P2 - m)} (17) 

is a non-trivial fact, particularly, because of the presence of the B-function in eq. (0). This 
G-function provide a correct kinematics in variable xn but cuts out a part of the integration 
domain in d^p. This cutting of the integration interval in the polar angle 6 leads to non- 
zero contribution of the matrix element containing 73, which is proportional to cos6. The 
sum rule (|T6|) for the first moment of is of a different nature, it presents the nucleon 

contribution into the total momentum of the deuteron [jlSj |T^ and ^at is a part of the total 



momentum carried by the non-nucleon component (mesons). An important pioperty of sum 
rules (|15D and ([l6| ) is that their r.h.s. do not depend upon the deuteron spin orientation. 
The SFs and 

6^2 are now calculated as follows: 




Mn/m 



F?M-')\_ [ .J/S) UfKe^). (19) 







where distributions f^l^ and A/^/'^ are given by 

f'^'iy) = lT.fTiy\ (20) 

Ar/^(l/) = /f/^(l/)-/o"^/^(l/). (21) 

The sum rules for f^^^{y) and Af^/^[y) follow from sum rules for fM^{y) ^'^^ defini- 
tions (0)-(|2l|): 

Mu/m 

I f"'{y)dy = \Y.{D\Q\D)M = i, (22) 
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Mo/m 



1 yr/^{y)dy=W{D\{&N)';\D)M = l-SN, (23) 



/m 

J Af''/''{y)dy={D\Q\D)M=i-{D\Q\D)M=o = 0, (24) 



Mo/m 

yAf"'{y)dy = {D\ {Q^f^ \D)m=i - {D\ (0^.)^ \D)m=o = 0- (25) 



Sum rules for the deuteron SFs 6f and 6^ are the immediate result of combining eqs. (p4|) - 
ID and dTD-dD: 



j dxDbfixD) =0, J dxDb^ixD) = 0, (26) 



i.e. in a full agreement with the sum rules (|1]) and (|^). 

An explicit expression for the distribution function fM^iu) {dJid therefore of f^^^{y) 
and Af^^^(y)) in terms of the components of the Bethe-Salpeter amplitude (eq. ([T^)) is 
quite cumbersome and it will be presented elsewhere. 



4 Nonrelativistic formulae 

The nonrelativistic expressions for f^^^{y) and Af^^^(y) can be obtained by using an anal- 
ogy of the charge densities calculated within the Bethe-Salpeter formalism and corresponding 
densities calculated with wave functions, i. e.: 

/ ^Tr{^'M(po,p)7o^Af(po,p)(p2-m)} ~ ^l,(p)^Af(p), (27) 

— oo 

[ ^Tr{^r,,(po,p)73*A/(Po,p)(P2-m)} ~ ^vp1,(p)v1>m(p), (28) 



2M, 



D 



-oo 



where ^m(p) is a nonrelativistic wave function of the deuteron (do not confuse with the 
Bethe-Salpeter amplitude, \E'Af (poP)') The well-known result for the spin- independent dis- 
tribution is immediately reproduced (see e.g. 0, O, 113): 



^ ^ M 



(27r)3 V m 



^^^co^ P (29) 



where u and w are the 5— and D-wave components of the deuteron wave function, respec- 
tively. The presence of the 0-function in the r.h.s. of eq. ( p9|) shghtly violates the sum rule 
(p2D. However, this usually phenomenologically is not noticeable, since the only region of 
large momenta, \p\ > 0.7 GeV, is affected by the 9-function and it does not contribute much 
to the norm of the deuteron wave function. 
For distribution Af^^^{y), one gets: 

Pl3 



m ^ J 
'''' s(y-^^^^^^^)Qiy) 



(27r)3 m 
(1 + \P^)P,(cose)^wi\p\) {2V2ui\p\) +w{\p\)}, (30) 

where P2 is the Legendre polynomial. 

Again, the sum rule ([2^ ) is broken by the presence of the 6-function in (^0]). Neglecting 
it, one gets 

1 Mo/m 



/ d{cose)il + ^^^)P2(cos^) = 0, (31) 
J m 



oc 



-1 



where the orthogonality property of the Legendre polynomials is used. 

A deviation from zero, caused by the 9-functions is not large compared to 1, but anything 
is large compare to 0! One can artificially adjust formula (|30D to fulfill this sum rule. For 
instance, small corrections to the normalization of the both terms with M = 1 and M = 
can be made to satisfy the sum rule in form (p^. However, the situation with the second 
sum rule, (^ and (^61), is more difficult and can not be fixed by any simple adjustments of 
the normalizations. Similar to eq. ([3T|), it can be written (neglecting the G-function!): 

1 Mn/m 


/\t)\cosO 
d(cos9)(MD -w+ \p\cose)(l + — )P2(cos9) ^ 0. (32) 
m 



7 



Thus, there is no reason for this sum rule to be satisfied with the nonrelativistic distribution 
function (pO]). 

In conclusion to this section, the nonrelativistic formulae, in principle, violates the sum 
rules for the SFs ftfg- However, one can still hope that it will be a small effect, one not 
noticeable in practice. The illustrative numerical calculations are given in the next Section, 
this helps one to understand a quantitative side of the effects. 



5 Numerical calculations 

In this Section results of the numerical calculations are presented. The SFs of the deuteron 
6^2(3^) &re calculated within both the relativistic and nonrelativistic approaches. 



The relativistic calculations are based on the formulae (0), (p^)-(|21|) and utilizes the 
realistic Bethe-Salpeter amplitude for the deuteron calculated in ref. [1^. This amplitude is 
essentially just a different presentation of the amplitude obtained by Zuilhof and Tjon [|T2| . 
An important feature of the calculations is a numerical "inverse Wick rotation" of the ampli- 
tude. This has been done by expanding a non-singular part of the matrix elements ([I^ ) into 
a series in po/m, up to the fourth order. Singularities in the nucleon propagators have been 
treated exactly. Note, the numerical approximation made can potentially cause a violation 
of the exact sum rules. 

The nonrelativistic calculations, eq. (pOD, uses the realistic wave function of the deuteron 
in the Bonn potential [|r^. Another ingredient of the calculations, the nucleon SFs Q"^), 
is taken from ref. ||2^ at = 10 GeV^. The results are neither very sensitive to the partic- 
ular choice of the parametrization for the nucleon SFs nor to the Q^- dependence of them. 

The distribution functions Af^^^(y) are calculated and the results are presented in Fig. 0. 
The relativistic (solid line) and nonrelativistic (dotted line) give very similar behavior of the 
distribution function. Indeed, it is difficult to distinguish between them, not speaking about 
making definite conclusions. The third line in the Fig. |^ is given for an illustration, and 
presents yAf^^^{y) for the relativistic calculations. The calculation of the sum rules is 
more representative. To understand the scale of effects, which are discussed below, it is 
customary to define auxiliary quantities: 

Mo /m Mb / m 

J Abs {a f'/^iy)) dye:, J Abs {y A f'/^iy)) dy c:, O.U. (33) 



The Bethe-Salpeter and nonrelativistic Bonn calculations give the same result in (|33D, with 
accuracy of ~ 5%. Thus, the effective distribution functions, Af^^^, are an order of mag- 
nitude smaller than the usual spin- independent distributions f^/^ normalized on 1. This is 
not a very important circumstance, but it works against accuracy in the numerical calcula- 
tions, since is a difference of two functions normalized on 1 (M = 1 and M = 0). 
Numerically, the sum rule ([^) is satisfied both in relativistic and nonrelativistic calculations 
with good accuracy, despite the approximate numerical "inverse Wick rotation" and the 
discussion after eq. (|30|). The correspondent integrals are ~ 5 ■ 10~^ and ~ 3 ■ 10~^ and 
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they should be compared to the estimate (0). The sum rule (^) may be used to improve 
distributions Af^^^ by making integrals for f^^^ and fo^^^ exactly the same. However, 
this does not lead to a significant variation of results for SFs, except x ^ for b^{x). 

The behavior of the b^{x) at x — > deserves to be considered more closely, especially 
for numerical calculations, since the nucleon function [x) can be divergent at small x. 
Unfortunately it is impossible to estimate bf{0) exactly for the realistic SFs F-l^ . However, 
a contribution of singularity can be evaluated. Indeed, let us assume a singular behavior as 
Ff ^ C/x, then eq. (|18D leads for small x to 

Mo/m r Mo/m x^ 

b^x^O) ~ I / Af^/^iy)dy = ^l J - J \ Af^^^ iy)dy 

X I J 

~ ^-CA/^/^(0), (34) 

X 

where Z = in exact relativistic formula, but it can be a small number in numerical 
calculations or in the nonrelativistic formalism. Thus, the limit of the deuteron SF b^[x) 
as X ^ is a constant, but one has to exercise great care in performing the numerical 
computations, since any error leads to a divergent behavior at small x. In this context, an 
adjustment of norms of the two terms in formulae (|2l|) and (|30|) has a meaning of subtraction 
of the numerical error from 6f at small x. 

The situation with the second sum rule @) is quite different. Numerically it is broken 
more significantly than the previous one. Correspondent integrals are ~ 1-10~^ and ~ 3-10^^ 
for relativistic and nonrelativistic calculations respectively, i.e. about 0.7% and 2% compare 
to (|33|). Therefore, numerical approximations slightly damage the relativistic formula. It 
is attributed to the numerical rotation to the Minkowski space. An adjustment of the 
normalization, as it has been discussed, slightly improves the accuracy (to 0.5%). On the 
contrary, the result for the nonrelativistic approach is stable with respect to any adjustments, 
since it is defined by the formulae (P^j). 

The SFs bi and 6^ are calculated within two approaches as well. The results are shown 
in Fig. ^ a) and b). The behavior of the functions in Fig. ^ a) suggests the validity of the 
sum rule (|1]). At the same time, the nonrelativistic calculation for 6^ in Fig. ^ b) (dotted 
line) obviously does not satisfy sum rule @). The main difference of the relativistic and 
nonrelativistic calculations is at small x, where these approaches give different signs for the 
SFs. To illustrate the effect of the presence of the G-function under integral in nonrelativistic 
formula (PP]), the calculations have been done as well with a restricted interval of integration 
over The condition |p| < 0.7 GeV corresponds to the "softer" deuteron wave function, 
but makes sum rule ( |5^ ) exact. Correspondent SFs are shown in Fig. ^ a) and b) (dashed 
line). The result of this "experiment" is that the effect of 0-function is not quantitatively 
significant. It also does not affect the principle conclusion about the second sum rule 
but makes the defect slightly smaller. This is understandable, since the sum rule breaking 
terms in (|32|) is oc \p\cos6. 
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6 Summary 

The structure functions (^) '^^ deuteron have been studied within the covariant ap- 
proach based on the Bethe-Salpeter formahsm for the deuteron. It is shown that the non- 
relativistic convolution model results in an incorrect behavior of these structures at small 
X and violates the exact sum rules. The importance of accurate relativistic calculations for 
6f*2(^) is demonstrated. 
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Figure captions 



Figure 1: Triangle diagram 



Figure 2: A/^/^(y). Curves: the BS - solid, the Bonn - dotted, the BS multiphed on extra 
factor y - dashed. 



Figure 3: &i^2(^)- Curves: the BS - solid, the Bonn - dotted, the Bonn with cut - dashed. 
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Fig. |I[ A. Umnikov, Relativistic calculation. 
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^ A. Umnikov, Relativistic calculation. 
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Fig. ^ A. Umnikov, Relativistic calculation. 
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